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The transport properties of an interacting one-dimensional quantum dot capacitively coupled 
to an atomic force microscope probe are investigated. The dot is described within a Luttinger 
liquid framework which captures both Friedel and Wigner oscillations. In the linear regime, we 
demonstrate that both the conductance peak position and height oscillate as the tip is scanned 
along the dot. A pronounced beating pattern in the conductance maximum is observed, connected 
to the oscillations of the electron density. Signatures of the effects induced by a Wigner molecule 
are clearly identified and their stability against the strength of Coulomb interactions are analyzed. 
While the oscillations of the peak position due to Wigner get enhanced at strong interactions, the 
peak height modulations are suppressed as interactions grow. Oscillations due to Friedel, on the 
other hand, are robust against interaction. 
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I. INTRODUCTION 

Quantum dot are an ideal playground to study the 
interplay between quantum confinement and Coulomb 
interactions. Beyond the well-known Coulomb blockade 
physics, lie several interesting physical effects which 
induce peculiar correlations among electron states. 
In two-dimensional (2D) quantum dots such correlations 
have been the subject of an intense theoretical research, 
especially employing several different numerical tech- 
niques.^ Of p articu lar interest is the emergence of a 
Wigner molecule the finite counterpart of a Wigner 
crystal.^ The high level of symmetry of typical 2D 
quantum dots (such as circular dots or pillars) results in 
rotating Wigner molecules, with a rotationally invariant 
density profile. As a result, the direct signatures of a 
Wigner molecule on the electronic density of the system 
are weak.lSIIIIIll Such states can be fully characterized 
only either considering their roto-vibrational spectrum^^ 
or by considering density-density correlation functions, 
whose experimental probe is particularly problematic. 
At the experimental level Wigner molecules have thus 
essentially been investigated by means of optical spec- 
troscopic techniques/^ while attempts at imaging the 
correlated electron wavefunction employing a scanning 
tunnel microscope (STM) have been put forward.^^ 
Due to the reduced dimensionality, interaction effects 
in one-dimensional (ID) quantum dots are even more 
dramatic. Indeed, interaction and quantum confinement 
effects are directly visible at the level of the electron 
density. Such ID quantum dots can be realized in 
several different ways, ranging from carbon nanotubes^ 
(CNTs) to cleaved edge overgrowth (CEO) wires. ^ 
The quantum confinement within a region of length L 
causes Friedel oscillations,^^ with a typical wavelength 
Xf = 27T/kF with kp the Fermi momentum. 
More intriguing is the formation of a Wigner molecule 
when Coulomb interactions exceed the kinetic energy. 



In ID it is pinned into the dot and induces a peculiar 
oscillatory pattern^^ on the electron density with a 
wavelength A^^ = Af/2. 

One-dimensional quantum dots have been investigated 
by means of numerical techniques, ranging from density 
functional approaches to exact diagonalizations,^^^*^ 
confirming the above picture. Despite their precision, 
such methods suffer of some limitation. They are usually 
restricted to a low number of particles and do not 
offer the flexibility of analytical models, which allow 
to investigate issues such as transport properties more 
easily. 

An analytical approach, widely employed to describe the 
low-energy sector of the physics of interacting ID elec- 
trons is the Luttinger liquid model. In connection 
with the bosonization technique, it represents a powerful 
method to deal with interacting ID fermions, allowing to 
explore the limit of not too low particle numbers. Inter- 
action effects are modeled by non-universal parameters 
gp and g^^ the strength of interaction in the charge and 
spin sectors respectively. This model has been applied 
extensively to the study of the transport properties 
of ID quantum dots.l^^lHSl j^^^g qI^q been applied to 
describe the formation of ID Wigner molecules. 43 46] 
Other models for ID fermions map onto a Luttinger 
liquid in their low-energy sector. The Hubbard m odel, 
for instance, can be mapped onto a Luttinger liquid,^^'^ 
with gp > 1/2, while g^ = I due to SU(2) symmetry. 
The extended Hubbard model removes the constraint 
on ^p.^^ Models for Wigner molecules consisting in anti- 
ferromagnetically coupled electrons oscillating around 
their equilibrium positions map^^ into a Luttinger 
liquid with = 1- 

The Luttinger model suffers some drawbacks. One is 
connected to the terms oscillating at wavelengths shorter 
than Xf in the series of harmonics^^ of its electron 
density. In general, the amplitudes of these terms 
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are model-dependent and are then free parameters!^ 
When including in the model terms describing Wigner 
oscillations, indeed all the amplitudes become weighted 
by phenomenological constants. Only by comparing 
with more refined methods one can attempt to determine 
such constants!^ 

In addition, treating both charge and spin in the Lut- 
tinger regime - the so called ^^spin-coherenf Luttinger 
model - is strictly valicP^^ only for temperatures and 
voltages smaller than the spin (and charge) bandwidth 
D(j = NiTVcr/L {Dp = NiTVp/L)^ where Vcr (vp) is the 
velocity of spin (charge) excitations. When such a 
constraint is not fulfilled, more refined models^^ should 
be employed including the so called ''spin-incoherenf 
liquid.i^^ 

Several methods are proposed and employed to ex- 
perimentally study Wigner molecules in ID. Besides 
spectroscopical tools, a Wigner molecule can be 
inferred from the modifications induced on the transport 
properties. Since Wigner oscillations are present in 
the electron density, it is also possible to directly probe 
such quantities in, e.g., momentum-res olved tunneling 
experiments with parallel quantum wiiesW^^^^^ 
Local probes are a promising technique. Recently, the 
injection from a STM tip has been theoretically proposed 
to detect local electron- vibron coupling,^ FriedeP^^'^ or 
even Wignei^^ oscillations. An STM is however sensitive 
to the tunneling density of states rather than to the 
electron density. More suitable is a charged atomic force 
microscope (AFM) tip, already proposed to image the 
spin-charge separation^^-^ in a Luttinger liquid. Due 
to the capacitive coupling to the electron density it 
allows to probe its oscillations. The effects of an AFM 
tip on the energy levels of a ID dot have been recently 
considered theoretically.^^ However, the infiuence on 
the conductance amplitude has not been addressed. 

In this work we fill the above gap, studying a ID quan- 
tum dot described as a Luttinger liquid, capacitively 
coupled to an AFM tip. The Luttinger model allows us 
to easily access the regime of large particle numbers, not 
yet considered.^^ We will consider the "spin-coherent" 
regime. In full generality, we will regard the interaction 
strength gp as a free parameter. The Friedel and Wigner 
oscillations of the electron density will be fully retained. 
We will develop a general and powerful framework which 
allows to systematically investigate both the linear and 
nonlinear transport properties to the lowest order in the 
tip interaction strength, in the requested regime. In this 
work, we will focus on the results concerning the linear 
transport regime. Our main findings are the following. 
Both the position and the height of the linear conduc- 
tance peak oscillate as a function of the tip position. 
While a shift of the position of the linear conductance 
peak has been already reported for small N^^^ the 
modulation of the height of the linear conductance peak 
is a novel result. These oscillations bring information 



about the Friedel or Wigner oscillations in the electron 
density. The oscillations induced by the Wigner molecule 
act differently on the conductance peak position and 
height as the interaction strength increases. In particu- 
lar, while the peak position oscillations due to Wigner 
get enhanced at strong interactions, the peak height 
modulations are suppressed. Oscillations due to Friedel, 
on the other hand, are robust against interaction both 
in the position and in the amplitude modulations. 

The scheme of the paper is the following. In Sec. |TI| we 
outline the model describing the system in the bosonized 



language. In Sec.|III]we evaluate the perturbative correc- 
tions to the dot chemical potential, while in Sec. [IV| we 



determine the tunneling rates in the presence of the AFM 
tip and employ them to evaluate the transport proper- 
ties. Analytic expressions for the linear conductance are 
provided in Sec. |lVl Section [V| contains our results . Fi- 
nally, in Appendix [A] we outline the evaluation of the 
tip-induced corrections to the tunneling rates. 



II. MODEL 

The system under investigation - see Fig. [l] - is an 
interacting one-dimensional (ID) quantum dot of length 
L capacitively coupled to a negatively charged {Vt < 0) 
atomic force microscope (AFM) tip and tunnel-coupled 
to source (S) and drain {D) contacts. A gate contact, 
biased at Vg and capacitively coupled to the dot, is also 
included. 

The system is described by the Luttinger model with 




FIG. 1. (Color online) Schematic setup of the ID quantum 
dot, perturbed by the negatively charged AFM tip (Vt < 0) 
at position xq , capacitively coupled to the gate contact poten- 
tial Vg and connected, via tunnel barriers at positions = 0, 
Xd = L with amplitudes ts,D, to source S and drain D con- 
tacts, kept at a potential difference V. 



a free band linearized around the Fermi point kp^s = 



Ns^^tt/ where Ng"^^ is the reference number of electrons 
with spin z component 5 = ± in the ground state. We 
consider here a reference state with an even total number 
= N^^ and introduce kp = kp^s- The odd case can 
be reproduced adding or subtracting one electron to the 
above case. 



3 



The dot Hamiltonian reads (from now on, h = 1) 



witlP 



n„>0 



(1) 

(2) 
(3) 



Here, Hjsj represents the contribution of the zero modes 



Np,„ = N. 



s=+ 



(4) 



with Ns the number of extra electrons with spin 5, with 
respect to Ns^\ The energies E^, = 7rVj^/2Lgi, have 
been introduced in terms of the velocity Vj^ of the mode 
v = p^a and of the Luttinger parameters g^. For re- 
pulsive interactions one has Qp = g < 1^ while g = 1 
corresponds to the noninter acting limit. On the other 
hand, g^j = 1 for an SU(2) invariant theory!^ The ve- 
locity of the charged mode is renormalized by the inter- 
actions and leads to Vp = vp/g where vp is the Fermi 
velocity. Within the standard Luttinger theory, the spin 
velocity is Vcr = vf even if, in more refined models, it 
may differ from this value j^^ ' ^^' Furthermore, Ng is the 
number of charges induced by a gate electrode capaci- 
tively coupled to the dot. 

Note that the parameter Ep might deviate from the sim- 
ple expression quoted above due to e.g. screening induced 
by surrounding gates. Therefore, Ep will be in the fol- 
lowing treated as a free parameter with Ep ^ E(j. 
The term describes collective, quantized charge and 
spin density waves with boson operators d^^nq and = 
TiVy/L^ with Uq G N. 

The electron field operator ^s(x) satisfying open bound- 
ary conditions ^s(O) = ^s(i^) = is 



(5) 



where ips,r{x) are 2L-periodic fermion fields represent- 
ing right (r = +) and left (r = — ) movers in the 
dot. Due to the open boundaries conditions one has 
'0s,r(^) = —i^s,-r{—x). The right movers operator ad- 
mits a bosonic representatiorF^ 



e'^^e V2 . (6) 



Here, a is the cutoff length, of order satisfies 



and r]s fulfill rjgTjs 



allowing the right anticommutation relations for different 
spins. The boson fields ^p(x), ^a{x) are given by 



-aq/2 



(cos {qx) - ig^ sin {qx)) dl , 



h.c. 



with q — UgTr/L. 



The coupling with the AFM tip is modeled as a capac- 
itive interaction 

Htip = -e dx Vtip{x)p{x) (7) 
Jo 

between the electron density —ep{x) and the tip potential 

Vtip{x)^ peaked around xq. 

The particle density p{x) = ^^^^ ps{x) with 

p,{x) = ^l{x)^s{x) (8) 
can be bosonized^^ using Eqns. (HJg), one has {a = k^^) 



Ps{x) 



(9) 



Here, the first three terms represent the long wave part 
expressed in terms of the antisymmetric field 

(Pp{x) ^s(p^{x) 



(Ps{x) 



V2 



while the Friedel part pf{x) is given 



I44T721 



with 



C{n, x) 
h{x) 



kp 

TT 

2kFX 



cos [C{Ns^x 
2i\x 



2^s{x)\ 
2h{x) . 



1 



tan 



sin(27rx/L) 



oTva/L 



■ cos(27rx/L)_ 



(10) 

(11) 

(12) 
(13) 

(14) 



In addition to the above terms, we will consider the so 
called Wigner contribution 

which arises from the presence of electron-electron in- 
teraction^^ beyond the Luttinger approximation, exter- 
nal perturbations,!^ or by effects of ion-electron interac- 
tions.*^ Its bosonic form iJ^ 



2C{Np/2,x)-2y/2ipp{x) 



(15) 



where A is a model dependent constant. Note that, in 
contrast to the Friedel term, the Wigner one depends on 
the charge sector only. The coupling with the tip is then 
represented as Hup = H^^ + H^^^ with 

H^^p = Vf^cos[C{Ns,xo) - 2^s{xo)] (16) 

s 

= Vw cos \2C{Np/2, x^) - 2V2^p{x^)\ , (17) 

where Vp and Vw are free parameters that depend on 
the shape of the AFM tip and on the weight of the 
Friedel/ Wigner oscillations. Note that we neglected the 
coupling with the long wavelength part of the density 
since it can be adsorbed in the Hamiltonian with a 
unitary transformation. 
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III. CHEMICAL POTENTIAL 

The presence of Hup shifts the energy levels of the 
quantum dot 



with respect to the bare case 
E°{Np,N, 



^{N,-N,f + ^Nl 



(18) 



(19) 



with (...) the thermal average with respect to at fixed 
particle number. The average is decomposed in Friedel 
and Wigner terms, for T ^ one has 



{Hup{xo))= Eip{Np,N,,xo) 

i=F,W 



(20) 



with 



EfipiNp, N^, xo)=VfK^{xo) cos [£(iV„ xo)] , (21) 



E^p{Np,N^,xo)=VwK'^{xo)cos[2C{Np/2,xo)] , (22) 



where K^{xo) = exp[—2{(pl{xo))] 

K^{xo) = exp[-4(^2(3.p)^] ^jth 



K^ixo) 



K'^ixo) 




n 2g 



sinh^ (if) +sin^ {^] 



and 



(23) 



(24) 



The chemical potential for a given configuration with 
Np charges and spin N^r is defined as 

l^d{Np, N„, xo) = E{Np + 1, iV^ ± 1, Xo) - E{Np, N„, xo) 

(25) 



with E{N) given in Eq. (18). The above expression holds 
either for ground states, where Np = {Np = 1) for even 
(odd) N {N being the total number of particles) and 



\Na\ = 



(-1) 



No 



(26) 



and for excited states where |TVcr| attains larger values 
with possibly Np > (TVp > 1) for an even (odd) N. 
The chemical potential can be decomposed as 



l^d{Np, TV,, Xo) = ^o{Np, TV,) + 6/i{Np, TV,, xq) , (27) 



where 



+ N,-N,]+E, 



l±2iV^ 



(28) 



lio{Np,N^) = Ep 



is the bare dot chemical potential and 



6fi{Np,N,,xo)= ^f^H^p.N^.^o) (29) 

are the corrections due to the tip. They have been the 
subject of numerical investigation^^ with exact diagonal- 
ization techniques in the regime of low TVp. As can be 
seen in Eq. (21) and Eq. (22), (5/i^(TVp, TV,, xq) exhibits 
an oscillatory shape enveloped by K^{xo). 
Let us now specify the above general expressions to the 
case involving ground states only, i.e. when Eq. (26) 
holds, relevant to study the linear transport regime. By 
exploiting the relation 2TV5 = TVp + sTV, it is easy to show 
that 



6/1^ {xo) = VfK^{xo) {cos [CiAf + 1, Xo)] - cos [C{J\f, xq)]} 



6fi'^{xo) = VwK'^{xo)<cos 



2jC 



TV. 



-,Xo 



cos 



2C{ ^,xo 



(30) 
(31) 



with 

if TV, is even 
\{Np -l)/2 ifTVpisodd ^ ^ 

The corrections to the chemical potential present an oscil- 
latory behavior given by the superposition of two cosine 
terms. For the Friedel case, if TV is even they oscillate 
with wavelength 2L/(TV + 2) and 2L/TV, while for odd 
TV their wavelength is 2L/(TV + 1) and 2L/(TV - 1). For 



the Wigner case, one always finds the spacial frequencies 
L/{N^1) and L/N . 

These oscillating patterns are modulated by the func- 
tions K^{xo) and K^{xo). They depend both on the 
number of electrons through the cut-off a = k^^ and on 
the interaction parameter ^p, dictating their power-law 
scaling. 

The interaction parameter g dictates the relevance of the 
Wigner term in comparison to the Friedel one. By in- 
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specting Eqns. ( 23|24 ) it is easy to show that Wigner 
correlations become relevant for ^ < 1/3, see Ref. 45, 
In addition, as N grows, a and thus K^{xo) are sup- 
pressed. In this limit quantum confinement effects be- 
come less relevant as the system crosses over towards the 
semi-classical regime. Furthermore, as the high-density 
limit is approached kinetic terms become more relevant 
than Coulomb repulsion, leading to a further suppres- 
sion of the signatures due to the Wigner molecule. 
The above analysis is confirmed by Figs. [2] and |3] which 
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FIG. 3. (Color online) Same as in Fig.[2]but for 16 ^ 17 and 
a — kp^ — 7t/8L. 



amplitudes of the tunnel barriers. The dot is subject 
to a symmetric voltage drop Vs = — V/2, Vd = V/2 
between the two leads. The leads are modeled as 
non-interacting Fermi gases with Hi^ads = Hs + Hjj 
and we have set the bare leads chemical potential to zero. 



FIG. 2. (Color online) Chemical potential corrections S/j^(xo) 
(units Ep) as a function of xq (units L) for the transition 
5^6 and different values of the interaction parameter: (a) 
g=l; (b) g=0.7; (c) g=0.4; (d) g=0.1. Red (gray) curves 
represent the Friedel correction, blue (dark gray) curves the 
Wigner correction. In all panels, Vf = Vw = and a = 
k-^ = L/Stv. 



show the corrections to the chemical potential for = 5 
and Np = 16 respectively, and different values of g. In- 
deed, for each of these cases it is clear that Wigner os- 
cillations grow and eventually become relevant as ^ ^ 0. 
Also, it is clear that for a given interaction strength both 
Friedel and Wigner oscillations get smaller as the number 
of partiles increases. 



IV. TRANSPORT 

The coupling between the dot and the leads is pro- 
duced by tunneling barriers at = and xjj = L 

Ht= E M,A(^A)V^.,+ (:rA) + /l.C.] =^i7^ 

X=S,Ds=± A 

with Xs,a(^a) = S^D) the fermion operators at 
the end point of the lead A and tx the transmission 



The key quantities of transport in the sequential 
regime are the tunneling rates. They will be evaluated 
via the time evolution of the density operator. 

At the initial time t = the dot is characterized by the 
state = \Np^Nl.). Here, we will consider the full dy- 
namics of the zero modes while the bosonic excitations 
will be assumed in thermal equilibrium. The initial den- 
sity operator then reads 



P(0) 



ZsZd ^6 



(33) 



with the partition functions = Tr {e ^^^}, r G 
{b,S,D} and P = l/{kBT). The probability Vs^^sf{^) 
of finding the dot in the final state \S^) = \N^^Nl) at 
time t is then obtained tracing out the leads and the 
bosonic degrees of freedom 



V. 



{t) = Tr{S^\p{t)\Sf). 



(34) 



with p{t) in the interaction picture^"^ with respect to 
p(t) = We-'^o'^^^^^A p(0) [fe^/o^^-V(-)l . (35) 



Here, V(r) = e'^ ^{Hup + Ht)e-'^ ^ and T (T) the 
time (anti-time) ordering operators. The probability in 
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Eq. (34) is computed to the lowest order in the tunnehng with V^i^^f{t) = 'P^[^_^^f{t) + ^V^i^^f{t) and 
barriers and in the couphng to the AFM tip. The typical 
structure is 



Vs^^sfit)= J2 (36) 

X=S,D I 



it)-- 



■- f dh f 
Jo Jo 



SV 



Jo 



{t)= / dh / dt2 / dt3J2F^^itut2,t3)+C■C. 



(37) 
(38) 



i=l,2 



where 

Fo\h,t2hTr{Sf\H,\h)p{0)H,\t2)\S^), 



(39) 



F,\h,t2,t3)=iTr{Sf\H^{h)p{0)H^{h)Hup{t2)\S^), (40) 
F^\h,t2,ts)=iTr{Sf\Ht\h)p{0)Hup{t3)Ht\t2)\Sf) . (41) 

Note that in the sequential tunneling limit considered in 
this paper the selection rules |Ar/ - 7V^| = |7V/ - 7V^| = 1 
apply. In the following, we will focus on the case where 
only two charge states are involved in the transport. 
The rates are given by the long time limit of the time 
derivate of V^{t) 



(42) 



Their explicit evaluation, within the bosonization formal- 
ism, is presented in detail in Appendix [Aj Here we quote 
the main results. The rates have the general expression 



^s^^sf — 
with 



nO,A 



ST 



F,A 



{xo)- 



-5tJ:^^,{xo). (43) 



representing the tunneling rate in the absence of the tip 
apart from the full chemical potential including tip cor- 
rections and 

(45) 

being the explicit corrections induced by the tip. Here, 
Tq = voltxf /na with the leads density of states and 
A^A = AneVx + fid{N;.K.xo) with An = (7V/ - 7V;) 
and iid{N' Ni.xo) defined in Eq. (p5|. The coeffi- 



cients Aq^^q^ and R^^\(^xo) are respectively defined in 



Eq. (|AT9| and Eq. ( [A2T| . 

The above rates fulfill the detailed balance relation 



(46) 



The framework developed here is very general an d allo ws 
to address both linear and nonlinear transport .1^^^ In 



I 

the rest of the paper, however, we will focus on the 
linear transport regime which already shows a rich and 
interesting physics. The nonlinear regime will be the 
subject of forthcoming investigations. 

In the linear regime only two charge values N and N 
will be considered, with the correponding ground state 
spins. Namely, 



= 0,iV, =0) 

:0,7V, = ±1) 



1,N^ = ±1) (TV even), 
1,N^ = 0) {N odd). 

(47) 



Note that in the absence of magnetic field, states with 
No- = ±1 are degenerate. 

The standard expression for the linear conductance, ex- 
pressed in terms of the rates involving the above states. 



436] 



G 



where F 



►AT 



(48) 



Np^N' ^ shorthand notation for F^i^^/, 

Dn = [3 + (-l)^+^] /2 is the degeneracy of the dot 
ground state with N electrons and 



The behavior of the conductance as a function of the ex- 
ternal parameters will be numerically investigated in de- 
tail in the next section. Here we discuss a useful analyti- 
cal approximation valid at low temperatures around the 
maximum of the conductance peak which is centered^^ at 
jj.(i{Np^ Nct^xq) = for T = 0. As shown in Appendix [A} 
in this regime the tunneling rates can be approximated 
as 



7\^o)/(/i) ; 



pA 



7\xo)/(-/i) 



where f{E) is the Fermi function, fi = iid{Np^N^^XQ) 
is th e dot chemical potential and ^'^{xq) is defined in 
Eq. (A29). Employing the above expressions and the 
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detailed balance in Eq. (46), one can rewrite the conduc- 
tance as 



with 



A(a;o) 



1 + (Z)iv/L'jv+i)e'3M 
7^(xo)7^(xo) 



(49) 



(50) 



Equation (49) describes a resonance peak located at 
ksT . fD 



log 



AT+l 



2 V Dn 
When expressed in terms of Ng^ the resonance is at 



E 



(51) 



(52) 



with 



AT * = AT + - + ^ 
^ 2 2E, 



log 



D 



D 



No 



(53) 



and S/i{xo) given in Eq. (29). Therefore, the position of 
the conductance oscillates due to the tip-induced spatial 



dci^o) = COS 

9s M = sin 



'27r(A/'+l) 
L 

27r(A/'+l) 
L 



fluctuations of the chemical potential. 
Besides the oscillations of the peak position, also the am- 
plitude of the conductance peak exhibits a modulation 
depending on the location of the tip: on resonance, the 
conductance evaluates to 



rA(xo), (54) 



and thus explicitly depends on xq via the ter m A (xq). 
Such modulation has not been reported so far. 1^2^ The 
term A{xo) can be expanded to within linear terms in 
(^r^''^. Assuming symmetric tunnel barriers with Fq = 
To, one has A{xo) ^ Aq + Af (xq) + Af'(xo) with 



(55) 



and 



Aiixo) = 2AoV^KHxo) J2 FHxo)gUxo) ■ (56) 



X — 2h{x) 
X — 2h{x) - 



Here, 

- 2kFX 
2kFX 



2itM 



sm 



L 

2iiN 



X — 2h{x) 



L 



X — 2h{x) 



(57) 
(58) 



with JV deflned in Eq. (32). The functions have the 
same functional form aHer the replacements J\f ^ N 
and h{x) 2h{x). The weighting functions Fc,s(xo) are 
given by 



RESULTS 



In this section we discuss the behavior of the linear 



conductance in Eq. (48). The full numerical results will 



fRxo) = {cos{kxo) + COS [k{l - xo)]} ,(59) 



b: 



FRxo) = ^ {sin(fexo) - sin [k{l - xo)]} , (60) 



be interpreted with the aid of the analytical expressions 
developed in Sec. [iVl 

We start considering a weakly interacting dot. Figure [4] 
shows the conductance for g = 0.9 and three transitions 
5 ^ 6, 12 ^ 13 and 20 ^ 21. For this calculation 
we have included only the corrections due to Friedel 
oscillations as one can expect that the Wigner term 
represents a vanishing perturbation due to Vw ^ in 
this limit .1^ The main panels show the conductance with 
a colorscale plot, as a function of the tip position xq and 
with m = {mi, m2,m3,m4 }, G N and rui > 0, Bl^ ^g- Below the density plot the maximum of the hnear 

conductance, Gm, is shown as a function of xq. Two 
main features are observed, namely 

{i) oscillations of its position, 
(ii) oscillations of its height. 



deflned in Eqns. ( |A24|A25[ ) , A = (mi + m2 ) + (^3 + 



7724), k — 7r(mi — 7722 + ma — 1714) /L and the summations 
are extended over the set of m 7^ {0, 0, 0, 0}. 
It can be shown that A^(xo) = A^(L — xq) and that 
A^(L/2) = 0, therefore the amplitude modulations are 
symmetric around the center of the dot and vanish there. 
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FIG. 4. (Color online) Color-map plot of the linear conduc- 
tance G (units e^To/ea) as a function of Ng and xq (units L) 
and plot of the conductance peak maximum Gm as a function 
of Xq (units L) for a transition N ^ N+1 with (a) = 5; (b) 
N = 12; (c) N = 20. The dashed lines superimposed to the 
color-maps represent the calculated chemical potential traces 
- see Eqns. ( |29|3Q|31t - the ones in the amplitude plots signal 
the value of the bare conductance peak obtained in the ab- 
sence of the tip. Here, only the Friedel corrections have been 
included, with Vf = 0.07ea and Vw = 0. Other parameters 
are ksT Omsa, Ep 5Ea and g 0.9. 



The fact (i) is in agreement with previous findings .1^11^ 
The fluctuations of the conductance peak position turn 
out to be proportional to the tip-induced correction 
to the chemical potential 5/i{xo) = 5/j.^{xo) in this 
case - see Eqns. ( 29|30 ). Indeed, the oscillations of the 
resonance position exhibit a number of maxima and 
minima in accordance with the discussion in Sec. IIIII 
This confirms the analytical prediction shown in Sec. |IV[ 
see Eq. (52). As a specific example, for the case of the 
transition 5^6 the plot in panel (a) exhibits three 
maxima and two minima. 

More interesting is the modulation of the height of the 
linear conductance which has never been reported so far. 
This modulation is sizeable and could be detected in a 
transport experiment. It exhibits a symmetric, beating- 
like pattern. The beating phenomenon is particularly ev- 
ident as increases. This fact will be analyzed in more 
detail later in this section. 

Let us now turn to the case of stronger interactions 
and study the effects due to the Wigner corrections. 
Figures [5] and |6] show the linear conductance for g = 0.3 
and the transitions 5^6 and 12 ^ 13 respectively. 
Panels (a) and (b) show the Wigner and the Friedel 
contributions to the conductance. While the Friedel 
corrections follow the pattern discussed above, Wigner 
corrections exhibit a shorter wavelength with precisely 
+ 1 maxima and minima, double the number of 
those observed in the Friedel case. The position of 
the conductance peak in the Wigner case follows the 
position-dependent correction to the chemical potential 
(5/i^(xo) defined in Eq. (31). Friedel and Wigner 



corrections induce therefore oscillatory patterns with 
very different wavelengths, allowing to distinguish the 
two mechanisms in a no n- ambiguous way. Although we 
have shown the two contributions separately, it can be 
in general expected that the Wigner and the Friedel 
mechanisms co-exist. Therefore, in panels (c) we show 
the conductance in the presence of both perturbations, 
considering the case of equal amplitudes. Even if here 
we have chosen the same weights for the Friedel and 
the Wigner terms (Vf = Vw)^ the latter impresses its 
peculiar oscillatory pattern with wavelength ^ L/N on 
the conductance peak position. Indeed, as we have al- 
ready observed in Sec. |III[ when interactions increase the 
Wigner contribution to the chemical potential becomes 
even more relevant than the Friedel one. This explains 
the persistence of A^ + 1 maxima and A^ minima in the 
oscillations of the chemical potential. Also conductance 
amplitudes exhibit a pattern reminiscent of the N-\-l and 
A^ maxima and minima but less pronounced with respect 
to the one shown in the position of the conductance peak. 

Let us now investigate the stability of these findings as a 
function of the strength of interactions. Figure [7] shows 
the linear conductance maximum and position for the 
transition 5 ^ 6 as a funciton of xq for several values 
of the Luttinger parameter g. As interactions grow, and 
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FIG. 5. (Color online) Color-map plot of the linear conduc- 
tance (units e^To/sa) for the transition A/" = 5 ^ 6 as a 
function of Ng and xq (units L) and plot of the conductance 
peak maximum Gm as a function of xq (units L) calculated 
with (a) Vf = 0, W = 0.076^; (b) Vf 0.07sa, Vw 0; 
(c) Vf — Vw — O.OTSfj- Other parameters are ksT — 0.05£cj, 
Ep and g 0.2. 



^ — ^ 0, two striking features are observed. While the 
oscillations of the chemical potential always exhibit six 
maxima and five minima, with an increase in the peak- 
to- valley ratio, the oscillations due to the Wigner contri- 
bution in the peak amplitude tend to vanish. 
Therefore, the effect of the Wigner corrections to the 
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FIG. 6. (Color online) Same plot as in Fig. [5] but for the 
transition AT = 12 ^ 13 for (a) Vf = 0, W = 0.07; (b) Vf 
0.07, Vw — 0; (c) Vf — Vw — 0.07ea- Other parameters are 
keT Omsa, Ep 5Ea and g 0.2. 



chemical potential and to the peak amplitude are very 
different: in the former case Wigner contributions are 
stable against Coulomb interactions, while Wigner cor- 
rections to the conductance maximum become vanishing 
for ^ ^ 0. 

In order to interpret the oscillations of Gm, let us turn to 
the analytic model developed in Sec. [TV| The height of 
the conductance peak follows A^(xo), given in Eq. (56). 
It is composed by the sum of two beating patterns, the 
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0.01 



-0.01 



Xq 




FIG. 7. (Color online) For the transition 5 ^ 6 (a) Conduc- 
tance maximum Gm (units e^To/so-) and (b) position offset 
Ng — Ng - see Eq. ( [52| - of the conductance peak as a function 
of Xq (units L) and several values of the interaction parameter 
g = 0.2, 0.1, 0.05, 0.02 (from darker to lighter color). Here 
we have chosen Vf = Vw — O.OTscr- Other parameters as in 
Fig-IU 



FIG. 9. (Color online) (a) Plot of (xo) (units e'^) 
as a function of xq (units L) for different values oi g = 
1, 0.8, 0.6, 0.4, 0.2 (from darker to lighter color). The arrow 
denotes the evolution of the family of curves as g decreases 
(interaction strength increases); (b) Same as in (a) but for 
{xq)] (c) Same as in (a) but for Ff (xo); (c) Same as in 
(a) but for F^ (xq). The dashed curve represents the asymp- 
totic case ^ ^ 0. All panels have been plotted for = 16. 



terms ^^^(xo), enveloped by the slow functions F^^^xq). 
The beatings arise from the superposition of the two co- 




FIG. 8. (Color online) (a) Plot of g^ (xq) as a function of xq 
(units L) for = 16 and a — kp^ — ii/^L] (b) Same as in 
(a) but for gY [xq). 



sine or sine terms occurring in ^|s(xo) with a wavelength 
difference of 27r/L and is therefore especially clear when 
N ^ 1. Figure [S] shows a typical example for the case 
of the Wigner correction. The Friedel case is perfectly 
analogous, only with half the wavelength. 

The envelope functions F^ ^{xq) parametrically depend 
on g and only very weakly on N through the cutoff a. 
They are plotted in Fig. [9) The overall shape of the 
Friedel and Wigner cases is analogous. Important dif- 
ferences however arise when the interaction strength in- 
creases {g 0). While F^^^xq) tend asymptotically to 



a limit function, the terms F^^(xo) 0. The math- 
ematical origin of this fact lies in the factor A~^ into 
the summations defining Eqns. ( 59|6Q ). Indeed, see Ap- 



pendix Al for the Wigner case A (x Sp (x g and thus 



Fj^(xoj~^ as ^ ^ 0. On the other hand, for the Friedel 
case A depends both on Sp and Sa- and thus non- vanishing 
contributions are still present even in the regime of very 
strong interactions. 

This confirms that the effects of the Friedel oscillations on 
the conductance peak amplitude seem more robust again 
the ones induced by Wigner oscillations. The rather 
counter-intuitive conclusion that the corrections to the 
maximum value of the linear conductance due to the 
Wigner term vanish in the ^ ^ limit is very inter- 
esting. We want to remind however that our discussion 
here concerns the intrinsic dependence on the interaction 
parameter, regardless the weighting factors of p^'^(x). 



VI. CONCLUSIONS 

In this paper we have calculated the transport prop- 
erties of an interacting one-dimensional quantum dot, 
described by means of the Luttinger liquid theory, in the 
presence of an AFM tip capacitively coupled to it. We 
have considered both Friedel and Wigner corrections to 
the electron density and evaluated their contributions 
both to the chemical potential and to the tunneling rates 
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to the lowest order in the AFM-dot couphng. 
To discuss the transport properties we have focused 
on the hnear regime, where we have shown that the 
AFM tip induces a shift of the conductance peak and 
a renormahzation of the conductance strength, which 
has never been discussed in previous works on this 
subject. Scanning the tip along the dot allows to observe 
oscillatory patterns related to both the Friedel and the 
Wigner oscillations of the electron density. A beating 
pattern emerges in the linear conductance height oscil- 
lations. The Friedel and the Wigner contributions show 
markedly distinct wavelengths allowing in principle to 
distinguish the two effects. Surprisingly, we have found 
that as the interaction strength grows, the effects of the 
Wigner oscillations on the conductance height vanish, 
while those on the oscillations of the peak position are 
reinforced. 

Interesting effects are expected to show up even in the 
nonlinear transport regime, also in connection to the 
possibility to trigger and probe dot excited states. This 
will be the subject of forthcoming investigations. 

Acknowledgments. Financial support by the EU- FP7 
via ITN-2008-234970 NANOCTM is gratefully acknowl- 
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with An = (TV/ - N'^) and fio{Nj^, N^) is the bare chem- 
ical potential of the dot, see Eq. (28). 



Note that we have introduced the full chemical potential, 
including tip-induced corrections. The kernel 



(A2) 



is the correlation function for lead A, which can be writ- 
ten as 



/oo 
-CO 



(A3) 



with ux the density of states in lead A and 

is the Fermi function. In the following we will assume 
v\ = uq. Concerning the dot, e"^^^'^^ stems from the 
thermal average over the collective excitati ons of th e ini- 
tial states and sum over the final ones witlPIMIlQl 



Wd{T) = 



1 



Wiep,T) + W{e,,T) 



(A5) 



Appendix A: Tunneling rate 

In this Appendix we outline the calculation of the tunnel- 



IV 



where \Sij) 



ing rates F^^^^j. introduced in Sec. 

|7v;./,ivy). 

Let us start with the term FQ^(ti,t2) in Eq. (39). Fol- 
low ing st andard procedures already outlined in litera- 
turd^^^' one obtains the standard result 



where 



where 



Wie,T) = J2 



-Tvan/ L 



n>0 

i sin(n£:r)] . 



(She 

coth ( ) (1 — cos{ner)) 



(A6) 



The perturbation terms 2^(ti, ^2, ta) in Eqns. (40 
can be decomposed into Friedel = F) and Wigner 
W) parts 

Ff{hMM)= Ff^{hMM). (A7) 



41) 



(Al) with 



■na 



^Ki{xo)F^{ti,h) cos CiiXx) + Agl'\tut2, h) 



(A8) 
(A9) 



where Xs = xq, Xd = L — xq, K^{xo) are defined in The functions 
Eqns. ( |23|24| ) and 



/^l^{x) = AnC{Nyf,x); 



i/f 



di^{x) = 2AnC{Nyf/2,x). 



(AlO) 
(All) 



iAgf\tut2,h) = {-iyg^{t3-h) 
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with j e {1,2}, are expressed in terms of 



where T{x) is the Euler gamma function. Inserting 
Eq. (IA8|) and Eq. (lAgl) into Eq. pTl and Eq. (|38|) it 



Q'^(t) = — - pW ( Ejy r + ) (A12) possible to perform the time integration exploiting the 

2 j^^p^p^^i \ J above Fourier expansions. 



The tunneling rate in Eq. (42) can be written as 



(A13) 



Since W{u) is periodic, it is convenient to exploit the 
Fourier series 



(A17) 



with 



V hf^^e-''^ . (A14) = To Ap,../ [A^A + ^P^P + (IcrSa] (A18) 



1 — — 00 



For /crT <C £(j one can approximate VF/ with their ex- , \, a , . 

pressioiPW calculated for T = 0, ^^""^^ = l*^' 



(1 



r(i + K)6i(0 
^ nr(i + K-o 



(A15) 



4 - U-,g/2i^-,l/2 . 



i\r{K) 



0{l), (A16) the corrections induced by the tip read 



where 



(A19) 



V^K^ixo) {cos [>C|(a;o)] - cos [4(a;o)] } (A20) 



(^n,m|A7^0 n,m 



1 ~ ^A,Q ^A,0 

A + 0+ ^ 



^(A)S. 



Fn(A^5,^p,^,),(A21) 
(A22) 



r 



with n = {ni, n2, ^3, 714}, m = {mi, m2, m3, 7714} and Furthermore, 



Cl!,{nH^) = COS \cf{x) + + fc^) 



A = Sp{ni + 712 + mi + ^2) + £a{ns + 714 + m3 + m4) , 
A = £p(mi + m2 - ni - 712) + ea{ms - ns - n^) , 

and kn = 7r(ni + 713 — n2 — n/^)/L (analogous for km)- 
(A23) Finally, the weig hts Bi ^ are 



^ ^-,l/25+,l/2^-,l/25+,l/2^+,l/25-,l/2^+,l/25-,l/2 
-^n,m "ni "^n2 "3 "4 nii "m^ "rm ' 



(A24) 
(A25) 



while ^ is expressed in terms of ^ as 

ni ,77,2 ,n3 ,724 , mi ,7712 ,7^13 ,7714 7T,i ,712 ,^3 ,714 ,7Tl2 ,''^1 ,''^^'4 ,'^3 
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The second term in Eq. (A20 



the first-order term in of t 



is easily recognized to be 
le expansion of 



where 



(A27) 



with iid{Np^Nl^XQ)^ containing the Friedel and Wigner 
corrections, defined in Eq. (25). From here o n, we will 



always include in both the bare rate, Eq. ( A26) and in the 



explicit corrections, Eq. ( A21 ), the full chemical potential 
including its corrections replacing AE^ with AEx. 
At low temperatures, j3£(j ^ 1, useful approximations to 
the tunneling rates when AEx ^ can be devised. This 
case is relevant for studying the linear conductance peak. 



Due to the argument of the Fermi function in Eq. ( A26 ) , 



the only term not exponentially suppressed is the one 
with Qp = Qa = 0. For the same reason, in Eq. (A21) the 



only terms which survive in the summations composing 
R^'^^q^{xo) are those with = = and ni = . . . = 
7^4 = 0. The tunneling rate can be thus written as 



where 



with 



and 



i=F,W 



r'o (1 



(A28) 



(A29) 



(A30) 



A 



-^0,mC'o,ii(^A) 



/(Ai?A) . 



(A31) 



^ L. p. Kouwenhoven, C. M. Marcus, P. L. McEuen, S. 
Tarucha, R. M. Westervelt, and N. S. Wingreen, in Elec- 
tron transport in quantum dots, NATO Advanced Studies 
Institute, Series E: Applied Science, edited by L. L. Sohn, 
L. P. Kouwenhoven, and G. Schon (Kluwer, Dordrecht, 
1997), p. 105. 

^ M. Dineykhan and R. G. Nazmitdinov, Phys. Rev. B 55, 
13707 (1997). 

^ P. Hawrylak and D. Pfannkuche, Phys. Rev. Lett. 70, 485 
(1993). 

^ M. B. Tavernier, E. Anisimovas, F. M. Peeters, B. Szafran, 
J. Adamowski, and S. Bednarek, Phys. Rev. B 68, 205305 
(2003). 

^ Y. Nishi, P. A. Maksym, D. G. Austing, T. Hatano, L. P. 
Kouwenhoven, H. Aoki, and S. Tarucha, Phys. Rev. B 74, 
033306 (2006). 

^ M. Rontani, C. Cavazzoni, D. Bellucci, and G. Goldoni, J. 

Chem. Phys. 124, 124102 (2006). 
^ A. Harju, H. Saarikoski, and E. Rasanen, Phys. Rev. Lett. 

96, 126805 (2006). 
^ C. Yannouleas and U. Landman, Phys. Rev. Lett. 85, 1726 

(2000). 

^ A. Puente, L. Serra, and R. G. Nazmitdinov, Phys. Rev. 
B 69, 125315 (2004). 

U. De Giovannini, F. Cavaliere, R. Cenni, M. Sassetti, and 

B. Kramer, New J. Phys. 9, 93 (2007). 

U. De Giovannini, F. Cavaliere, R. Cenni, M. Sassetti, and 

B. Kramer, Phys. Rev. B 77, 035325 (2008). 

F. Cavaliere, U. De Giovannini, M. Sassetti and B. Kramer, 

New J. Phys. 11, 123004 (2009). 

S. M. Reimann and M. Manninen, Rev. Mod. Phys. 74, 
1283 (2002). 



^ C. Yannouleas and U. Landman, Rep. Prog. Phys. 70, 2067 
(2007). 

E. Wigner, Phys. Rev. 46, 1002 (1934). 
S. Kalliakos, M. Rontani, V. Pellegrini, C. P. Garcia, A. 
Pinczuk, G. Goldoni, E. Molinari, L. N. Pfeiffer, and K. 
W. West, Nat. Phys. 4, 467 (2008). 

M. Rontani, E. Molinari, G. Maruccio, M. Janson, A. 
Schramm, C. Meyer, T. Matsui, C. Heyn, W. Hansen, and 
R. Wiesendanger, J. Appl. Phys. 101, 081714 (2007). 
M. Bockrath, D. H. Cobden, P. L. McEuen, N. G. Chopra, 
A. Zettl, A. Thess, and R. E. Smalley, Science 275, 1922 
(1997). 

A. Yacoby H. L. Stormer, N. S. Wingreen, L. N. Pfeiffer, 
K. W. Baldwin, and K. W. West, Phys. Rev. Lett. 77, 
4612 (1996). 

G. F. Giuliani and G. Vignale, Quantum Theory of the 
Electron Liquid (Cambridge University Press, Cambridge, 
2005). 

Y. Gindikin and V. A. Sablikov, Phys. Rev. B 76, 045122 
(2007). 

J. S. Meyer and K. A. Matveev, J. Phys: Condens. Matter 
21, 023203 (2009). 
2^ W. Hausler and B. Kramer, Phys. Rev. B 47, 16353 (1993). 
D. Agosti, F. Pederiva, E. Lipparini, and K. Takayanagi, 
Phys. Rev. B 57, 14869 (1998). 

G. Bediirftig, B. Brendel, H. Frahm, and R. M. Noack, 
Phys. Rev. B 58, 10225 (1998). 

B. Szafran, F. M. Peeters, S. Bednarek, T. Chwiej, and J. 
Adamowski, Phys. Rev. B 70, 035401 (2004). 

2^ E. J. Mueller, Phys. Rev. B 72, 075322 (2005). 

S. H. Abedinpour, M. Polini, G. Xianlong, and M. P. Tosi, 
Phys. Rev. A 75, 015602 (2007). 



14 



L. Shulenburger, M. Casula, G. Senatore, and R. M. Mar- 
tin, Phys. Rev. B 78, 165303 (2008). 
^° A. Secchi and M. Rontani, Phys. Rev. B 80, 041404(R) 
(2009). 

J. Qian, B. I. Halperin, and E. J. Heller, Phys. Rev. B 81, 
125323 (2010). 

G. E. Astrakharchik and M. D. Girardeau, Phys. Rev. B 
83, 153303 (2011). 

T. Giamarchi, Quantum Physics in One Dimension, Ox- 
ford Science Publications (2004). 
J. Voit, Rep. Prog. Phys. 58, 977 (1995). 
R. Egger and H. Grabert, Phys. Rev. Lett. 79, 3463 (1997). 
A. Braggio, M. Sassetti, and B. Kramer Phys. Rev. Lett. 
87, 146802 (2001); M. Carrega, D. Ferraro, A. Braggio, N. 
Magnoh, and M. Sassetti New J. Phys. 14, 023017 (2012). 
T. Kleimann, F. Cavaliere, M. Sassetti, and B. Kramer, 
Phys. Rev. B 66, 165311 (2002). 

J. U. Kim, L V. Krive, and J. M. Kinaret, Phys. Rev. Lett. 
90, 176401 (2003). 

F. Cavaliere, A. Braggio, J. T. Stockburger, M. Sassetti, 
and B. Kramer, Phys. Rev. Lett. 93, 036803 (2004). 

F. Cavaliere, A. Braggio, M. Sassetti, and B. Kramer, 
Phys. Rev. B 70, 125323 (2004); E. Paladino, A. D'Arrigo, 
A. Mastellone, G. Falci, New J. Phys. 13, 093037 (2011). 
L. Mayrhofer, and M. Grifoni, Eur. Phys. J. B 56, 107 
(2007). 

C. D. Graf, G. Weick, and E. Mariani, Europhys. Lett. 89, 
40005 (2010); G. Cuniberti, M. Sassetti, and B. Kramer, 
Phys. Rev. B 57, 1515 (1998). 

H. J. Schulz, Phys. Rev. Lett. 64, 2831 (1990). 
H. J. Schulz, Phys. Rev. Lett. 71, 1864 (1993). 

L Safi and H. J. Schulz, Phys. Rev. B 59, 3040 (1999). 
S. A. Soffing, M. Bortz, L Schneider, A. Struck, M. Fleis- 
chhauer, and S. Eggert, Phys. Rev. B 79, 195114 (2009). 

G. A. Fiete, K. Le Hur, and L. Balents, Phys. Rev. B 73, 
165104 (2006). 

K. A. Matveev, A. Furusaki, and L. L Glazman, Phys. Rev. 
Lett. 98, 096403 (2007). 

K. A. Matveev, A. Furusaki, and L. L Glazman, Phys. Rev. 
B 76, 155440 (2007). 

G. A. Fiete, Rev. Mod. Phys. 79, 801 (2007). 

K. A. Matveev, Phys. Rev. Lett. 92, 106801 (2004). 

F. D. M. Haldane, Phys. Rev. Lett. 47, 1840 (1981). 

K. Hiraki and K. Kanoda, Phys. Rev. Lett. 80, 4737 

(1998). 

P. Horsch, M. Sofin, M. Mayr, and M. Jansen, Phys. Rev. 
Lett. 94, 076403 (2005). 

L. H. Kristinsdottir, J. C. Cremon, H. A. Nilsson, H.Q. Xu, 
L. Samuelson, H. Linke, A. Wacker, and S. M. Reimann, 
Phys. Rev. B 83, 041101(R) (2011). 



V. V. Deshpande and M. Bockrath, Nat. Phys. 4, 314 

(2008) . 

^'^ H. Steinberg, O. M. Auslaender, A. Yacoby, J. Qian, G. 
A. Fiete, Y. Tserkovnyak, B. L Halperin, K. W. Baldwin, 
L. N. Pfeiffer, and K. W. West, Phys. Rev. B 73, 113307 
(2006). 

O. M. Auslaender, H. Steinberg, A. Yacoby, Y. 
Tserkovnyak, B. L Halperin, K. W. Baldwin, L. N. Pfeiffer, 
and K. W. West, Science 308, 88 (2005). 

G. A. Fiete, J. Qian, Y. Tserkovnyak, and B. L Halperin, 
Phys. Rev. B 72, 045315 (2005). 

N. Traverso Ziani, G. Piovano, F. Cavaliere, and M. Sas- 
setti, Phys. Rev. B 84, 155423 (2011). 
L. C. Venema, J. W. G. Wildoer, J. W. Janssen, S. J. Tans, 

H. L. J. Temminck Tuinstra, L. P. Kouwenhoven, and C. 
Dekker, Science, 283, 52 (1999). 

^2 S. Eggert, Phys. Rev. Lett. 84, 4413 (2000). 

A. Crepieux, R. Guyon, P. Devillard, and T. Martin, Phys. 
Rev. B 67, 205408 (2003). 

G. Buchs, D. Bercioux, P. Ruffieux, P. Groning, H. 
Grabert, and O. Groning, Phys. Rev. Lett. 102, 245505 

(2009) . 

S. Pugnetti, F. Dolcini, D. Bercioux, and H. Grabert, Phys. 
Rev. B 79, 035121 (2009). 

A. Nocera, C. A. Perroni, V. Marigliano Ramaglia, and V. 

Cataudella, Phys. Rev. B 86, 035420 (2012). 

A. Secchi and M. Rontani, Phys. Rev. B 85, 121410 (2012). 

J. Lee, S. Eggert, H. Kim, S.-J. Kahng, H. Shinohara, and 

Y. Kuk, Phys. Rev. Lett. 93, 166403 (2004). 

L. M. Zhang and M. M. Fogler, Nano Lett. 6, 2206 (2006). 

L Ussishkin and L. L Glazman, Phys. Rev. Lett 93, 196403 

(2004). 

E. E. Boyd and R. M.Westervelt, Phys. Rev. B 84, 205308 
(2011). 

'^^ M. Fabrizio and A. O. Gogolin, Phys. Rev. B 51, 17827 
(1995). 

Other derivations of the density operator lead to a form 
which only slightly differs at the dot boundaries from the 
one we use.^^^ 

K. Blum, Density Matrix Theory and Applications, Plenum 
Press (1981); M. Grifoni, M. Sassetti, and U. Weiss, Phys. 
Rev. E 53, R2033 (1996); E. Paladino, M. Sassetti, G. 
Falci, and U. Weiss, Phys. Rev. B 77, 041303 (2008). 

F. Haupt, F. Cavaliere, R. Fazio, and M. Sassetti, Phys. 
Rev. B 74, 205328 (2006); M. Merlo, F. Haupt, F. Cava- 
liere, and M. Sassetti New J. Phys. 10, 023008 (2008). 

G. Piovano, F. Cavaliere, E. Paladino, and M. Sassetti, 
Phys. Rev. B 83, 245311 (2011); M. Carrega, D. Ferraro, 
A. Braggio, N. Magnoli, and M. Sassetti, Phys. Rev. Lett. 
107, 146404 (2011); D. Ferraro, A. Braggio, N. Magnoh, 
and M. Sassetti Phys. Rev. B 82, 085323 (2010). 

'^'^ C. W. J. Beenakker, Phys. Rev. B 44, 1646 (1991). 



